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In this talk, I will discuss the predictions of the first results from RHIC: the charged 
particle multiplicity dN c h/dr}, its centrality dependence and the elliptic flow t> 2 . 



1. Charged particle multiplicity in central collisions 
1.1. Predictions 

To predict particle multiplicities in pp and AA collisions from first principles has so far 
been impossible. Consequently, various models have been introduced. In Fig. 1 I have 
collected predictions for the charged particle multiplicity dN^/dy at y — for Au-Au 
collisions at y/s = 200 AGeV, presented at Quark Matter '99 (]J and after, until July 
2000, i.e. prior to the release of the very first data from RHIC. From the top, the figure 
shows predictions from event generator models Heavy Ion Jet Interaction event Genera- 
tor 0, H1JWG+ Zhang's Parton Cascaded A Relativistic Transport model ||, Relativis- 
tic Quantum Molecular Dynamics || , Ultrarelativistic Quantum Molecular Dynamics |5] , 
VNI+UrQMD g], Hadron String Dynamics and VNI+HSD |7j, NEXUS @, Dual Parton 
Model ||, DPMJET ||10|| , String Fusion Model pTTH , and from other models, such as Linear 



Extrapolation of Ultrarelativistic nucleon-nucleon Scattering to Nucleus-nucleus collisions 



12] , EKRT saturation model Hydrodynamics+UrQMD [ Hfl , thermal fireball model 
i~5f and McLerran-Venugopalan model [[RJ. The assumptions of the underlying QCD dy- 
namics, and even the degrees of freedom (partons, strings, hadrons, classical fields) vary 
between different models, resulting in a variation of a factor two in the predictions. The 
multiplicities at yfs = 200 AGeV are yet to be measured but, based on the PHOBOS 
data at 56 and 130 AGeV [17 and an extrapolation to 200 AGeV (the region between 



the dotted lines), it seems that most models predict too many particles per unit rapidity. 

1.2. Comparison with the first data 

Given the large theoretical uncertainties, the appearance of the first RHIC data on 
dN c h/dr] from Au-Au collisions at y/s = 130 AGeV []T7J is greeted with a great enthusiasm. 



The model predictions can now be compared with actual measurements for the first time 
at collider energies. In Fig. 2, I have shown examples of such comparisons. In LEXUS 



18| , one first fits all the parameters at nucleon-nucleon level. Extrapolation to ultra- 
relativistic heavy ion collisions is then done by assuming different dynamics. As nicely 
illustrated in the figures, the "linear" extrapolation to AA, where every binary collision 
contributes to particle production, gives the best fit at the SPS, whereas the wounded 
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nucleoli model fits best at both RHIC energies. A suppression mechanism of the produced 
quanta is also needed in string-based models, such as in the Dual String Model |19| , where 
fusion of strings is needed to fit the RHIC data. HIJING p5| combines pQCD parton 
dynamics with that of strings. On one hand, a strong gluon shadowing in the hard 
pQCD component of HIJING reduces the multiplicities, but on the other hand this is 
partly compensated by jet quenching which feeds some of the jet energy into the system, 
resulting in the multiplicities close to the measured values ||20|| . The EKRT saturation 
model, which is based on pQCD parton dynamics supplemented by a requirement of 
saturation of produced gluons also predicted correctly both the normalization and the 
y/s scaling for central collisions, as shown by the dotted-dashed line in Fig. 2 (lower 
left). Another model where the observed multiplicities can be well explained is AMPT, A 



MultiPhase Transport Model pT[| , which takes the initial conditions for the partonic and 
hadronic afterburners from HIJING and fits some of the model parameters from the SPS 
data. Of the above models only ]13[ is a published prediction for 56 and 130 AGeV but 
to my knowledge also [fL^] , p0| (see Fig. 3) and |?T) are predictions in the sense that the 
model parameters were not tuned to fit the first RHIC data. 
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rigure 1. Predictions for dN^/dy 
in Au+Au at y=0, y/s = 200 ^GeV 
from different models before the ap- 
pearance of the first RHIC data. 
The information in the parentheses 
specifies the parameters used. The 
time of releasing the result is indi- 
cated on the left. An approximate 
factor 2/3 has been applied to con- 
vert iVtot (asterisks) to N<± (boxes). 
Factors 1.1 and 0.9 approximately 
accounting for the conversion of n 
into y and for a ~ 5% centrality se- 
lection, respectively, have not been 
applied in the figure. The vertical 
dotted lines are obtained by scaling 
the PHOBOS data at 130 ,4GeV for 
dN c h/dn (6% central) up by 1.1 * 
(200/130) - 37 based on the observed 
-^/s-scaling and to account for the 
conversion of n into y. In prepara- 
tion of this figure I have partly used 
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the review [11|. 



An implication from the comparison in Fig. 2 is that coherence phenomena in particle 
production become important at RHIC energies: particle multiplicities are less than what 
could be expected based on a mere linear extrapolation from SPS to RHIC based on pp 
physics. It also seems that efficient final state interactions are needed, which in turn 
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points towards an early formation of pressure and thermalization. It is also obvious that 
even within each model the theoretical uncertainties related to the fact that the underly- 
ing dynamics is not precisely known in advance, are clearly larger than the experimental 
(systematic) errorbars of the data. Therefore, the measured multiplicities are indispens- 
able constraints for the models. Essential further constraints, hopefully stringent enough 
to rule out some models, will be obtained from the other global observable, dEx/dr], which 
is more sensitive to the actual evolution of the system than dN^/dr]. 
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Figure 2. Upper left: LEXUS |18| (solid), wounded nucleon model (dashed); upper right: Dual 
String Model with (dashed) and without (dotted-dashed) string fusion |l~9[| ; lower left: HIJING 
with shadowing and with and without jet quenching (Au+Au, solid) f20[| t ogether with the 
EKRT saturation model (dotted-dashed) |]|; Lower right: AMPT modelpl|. 



1.3. Perturbative QCD in dN ch /dr] 

Towards higher cms-energies, perturbative production of gluons and quarks in AA col 
lisions becomes more important, eventually even dominant, as first suggested in P^-EB 



The agreement of HIJING and EKRT-type models with the first RHIC data suggest that 
at RHIC we are witnessing an onset of pQCD particle production. 

In first approximation (leading twist), perturbative production of partons, minijets, 
can be computed assuming factorization of the long- and short-distance physics, i.e. the 
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nuclear parton distributions f^(x,Q) and the partonic cross sections fry. Schematically, 
the integrated minijet cross section for producing partons with transverse momenta larger 



than a cut-off scale p 



1..2 GeV, is afj 



(Pt >Po)=K Eij ft ® if ® Oiy So far this 



has been computed in lowest order pQCD only and a iGfactor is needed to simulate the 
NLO contributions. The good news is, however, that due to the recent progress |27]], the 
iG-factor can now be computed exactly in certain cases, such as in Et production. The 
order of magnitude is K ~ 2, but the exact value depends on y'i and on the parton 
distributions chosen. 
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Figure 3. Top left: The HIJING prediction for dN^/dr] in Au-Au collisions at ^/s = 130 
j4GeV [29]. Bottom left: Determination of the dynamical saturation scale p sa t = po(A,y/s) in 
the EKRT saturation model [jig] . Ni is the number of produced pQCD partons in |y| < 0.5. 
Right: Comparison of the EKS98 nuclear effects (filled circles) [EH] and the HIJING default 
parameterization (open diamonds) with the ratio of structure functions F^/F^ 11 measured by 
NMC |3(| (open boxes with errorbars). 



The largest uncertainty in minijet production is the strong dependence on the minimum 
transverse momentum scale po allowed. In two-component models po determines the 
division into the pQCD and non-perturbative particle production. In e.g. HIJING [pql , 
Po = 2 GeV (with K = 2) is determined from the pp multiplicity data, and it is kept 
fixed in the extrapolation to AA. This makes the perturbative component to scale as 
N ~ A 4 / 3 . In the pQCD models supplemented by the requirement of gluon saturation 
0,0, po is determined dynamically as the scale which dominates all particle production 
in A A collisions. This leads to p sat = po(y/s,A), as shown in Fig. 3 (bottom left, the 
intersections of the curves) for the EKRT saturation model. 

For the modifications of the parton distributions in nuclei, ft(x, Q) 



R?(x,Q)fi(x,Q), 
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different parametrizations are available. The minimum requirement is that a relevant 
parametrization is consistent with the existing extensive data in deeply inelastic IA scat- 
tering. If collinear factorization is used, as usually is the case in models with minijets, the 
nuclear modifications should also be compatible with pQCD scale evolution. The EKS98 
nuclear parton distributions |[28|| , now included also in the CERN-PDFLIB, fulfill these 



requirements. An example of the constraints imposed by the DIS data is shown in Fig. |3|. 
The effects caused by different nuclear parton distributions in the computation of the 
multiplicities in HIJING is shown by Fig. ||] (top left). 

1.4. Gluon saturation in the initial and final state 

In modeling initial particle production in AA collisions at collider energies, more em- 
phasis has been recently given to models with gluon saturation. The idea can be sketched 
briefly as follows: let N g (Q, AY) gluons appear at a scale Q within a rapidity correla- 
tion length AY. Saturation of the number of gluons in the transverse plane takes place 
when N g (Q, AY) x a g {Q) ~ ttR\, where cr g (Q) ~ Oi s {Q 2 )/Q 2 is a cross section for 2 — > 1 
process and kR 2 a is the available total transverse area. Then, provided that N g (Q,AY) 
and (Tg{Q) can be computed, the saturation scale Q saXl = Q which fulfills the saturation 
criterion, can be found. In the case of heavy nuclei at high cms-energies, N g becomes 
large and saturation scales become Q sat ^> Aq CD . 

Gluon saturation in the initial state, in the wave functions of the colliding objects, 



was first discussed in |31] for pp and in |22| , |23| for AA collisions. The idea of initial 
state saturation is present also in the classical field approach P^J53], where Q 2 at = 



j^^ a s xG(x, Qg at )T^(g) with the gluon distribution xG and the nuclear overlap func- 
tion Ta(s). The saturation in the final state is obtained through the initial state gluon 
saturation, as the transverse profile of the rapidity density of produced gluons in a cen- 
tral AA collision can be expressed as |HJ dN AA /d 2 sdy = cCFQ 2 at /(a s 2n 2 ). The "gluon 
liberation constant" c is expected to be of the order of one P3| |. Analytical calculations 



in the semiclassical approximation |34| give c = 1.39, and it is interesting to note that if 
an agreement with the PHOBOS data is required, one gets Q 2 at = 2.1 GeV 2 Aq CD . 
The SU(2) lattice formulation of the classical field approach [[£J gives c = 1.29 ± 0.09, 
the SU(3) calculations are in progress. Similar values are also obtained in |36|], based on 
fits to the PHOBOS data. 

Saturation in the final state, i.e. of produced gluons is also possible, even without the 
requirement of saturation in the initial state. In the EKRT saturation model [|13] the 



number of produced gluons above a minimum p^-scale po, Naa{po, y/s, A), is computed 
from pQCD. The cross section in the saturation criterion is simply taken as the geometrical 
one, <j g = 7r/pg, not specifying any powers of a s or group theoretical factors. [] The 
saturation criterion in a central rapidity unit of central collisions, Naa{Po, \/s) x n/p^ = 
tcR\, then determines the dominant scale p sat = po(^/s,A) as shown in Fig. ||. At a 
scaling limit <Xj e t ~ Po 2 , the saturation scale behaves as p 2 &t ~ A 1 / 3 , causing the initial 
state multiplicity to grow as Naa ~ A instead of A 4 / 3 . A full calculation with realistic 
parton distributions accounting for the small-x increase and the EKS98 nuclear effects 
gives a multiplicity of produced partons as |T3] A^Gosat) = 1.383A°- 922 ( v /i) - 383 . 



lr The suppression of the production of small-pT quanta, leading to the saturation, can also be viewed in 
terms of a screening mass generated into the produced system fj9f . 
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Also the released transverse energy, E^f(p sat ) at saturation can be computed from 
pQCD. As the formation time of the produced system is r, ~ l/Ps&t, the initial densities 
in the EKRT approach can be computed through converting E^f(p SSbt ) into a thermal 
energy density e$, which further converts into a thermal number density nf 1 . On the other 
hand, conversion of the initial number of gluons at saturation, NAA(Ps&t), into rii gives 
densities very close to the thermal n\ h . In this sense the system looks thermal already at 
saturation [13|. Thus assuming thermalization at Ti, and an isentropic expansion stage 
described in terms of (boost invariant) hydrodynamics, the final state multiplicity of 
charged particles is computable from the entropy produced initially at 7$. This leads to 
a simple scaling law |]J dN ch /dy ~ § x ^fN t = 0.83 J 4°- 922 (y / i) - 383 for central collisions. 
Agreement with the PHOBOS data is amazingly good, as seen in Fig. 2 (lower left). More 
detailed description of the expansion stage in terms of hydrodynamics with transverse flow 
effects and resonance decays has now also been done |37| , the results are shown in Fig. 4. 



Figure 4. dN^/drj and dEx/dr] aver- 
aged in \rj\ < 1 vs. -y/s, as obtained 
from pQCD+saturation+2d hydrody- 
namics+resonance decays [37]. The 
open symbols are for central collisions, 
the closed ones are with a 6% centrality 
selection. In the left panel, the lines are 



from [13] for central collisions, the open 
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2. Centrality dependence of dN ch /dr] 



The dependence of dN^/drj on the centrality of the collision has been suggested in [20 



for testing the predictions from different models. Before the release of the PHENIX data 
on the centrality dependence, three concrete predictions were submitted: 

1. In HIJING p0| the multiplicity is formed as a combination of a component which 
scales as the number of participants ~ A, and a pQCD component which scales with 
the number of binary collisions ~ ^4 4 / 3 (keeping p fixed): dN ch (h)/dr] = Ap art (b)n soft + 
fN hm (h)a^( Po ) 



2. The EKRT saturation model [pl| can be extended to describe the saturation of gluons 
locally at a transverse distance s in a collision with an impact parameter b. Assuming 



the extreme case that saturation fully dominates the particle production, one obtains p8 
dN ch (b)/dr] w 0.9 x 0.9 x § x / d 2 sp 2 sat (^, A, b, s)/tt. 

3. Dividing the particle production into a "hard" and a "soft" components which scale 
A r bin, correspondingly, one may write |36| dN c ^(b) / dr] = xn pp N piLr t(b) + 

is obtained from the pp data, 

In H 



as 



Ap art and 



(1 — x)n P pAbin(b)/2. The average particle multiplicity in n, 
and the fraction x ~ 0.05 . . . 0.09 is obtained by fits to the PHOBOS data [17 
also an extension of the initial state saturation model 



to noncentral collisions was 
suggested in the form dN ch /dr] = ^cN p3irt xG(x, Q S at(b) 2 ), with c obtained from a fit to 
PHOBOS data. 

As shown in Fig. 5 (lower right), the centrality data measured by PHENIX ||D| lie 
in between the EKRT type saturation model and HIJING, quite close to the behaviour 
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Figure 5. Multiplicity dN ch /dr]/(0.5N pavt ) vs. iY part at 



y/s = 130 ,4GeV. From left: HIJING @ (histograms); lo- 
cal saturation model |38| (solid) and EKRT saturation |13| 
(dotted); "hard" +" soft" eikonal approach (crosses) and ini- 
tial state saturation model (diamonds) from |36fl . Lower 
right: the PHENIX data @. 
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obtained in . Obviously, gluon saturation (EKRT) is not the dominant mechanism for 
non-central (N pa _ rt ^ 200) collisions at \fs = 130 ^4GeV but a non-saturated component is 
needed at small transverse densities. For central collisions, however, gluon saturation is 
not ruled out as the dominant mechanism, as also indicated by the PHOBOS data. 

The multiplicity measurements alone are not sufficient to rule out models, as most mod- 
els can be tuned to fit the measured multiplicities and even their centrality dependence. 
For instance, HIJING has theoretical uncertainties in the division into soft and hard 
physics and in the extrapolation to AA, and also in the pQCD component as discussed 
above. The saturation models in turn suffer from uncertainties related to the saturation 
criterion which may be of non-perturbative nature and contain unknown constants and 
powers of (running) a s which may need to be phenomenologically determined. In any 
case, it is unlikely that the true theoretical errors in any of the models could be squeezed 
below the ±10% level already reached by the first RHIC data. 

Of the global variables, Et should be more efficient in ruling out models than iV c h, as 



Et is much more sensitive to the evolution of the system: as shown in |[13|j37| , the ratio 
of the initially released Et (at saturation) to the measurable final state Et is roughly a 
factor three. The reduction is due to the pdV work during the expansion stage. Such 
reduction does not take place in models which do not describe the final state interactions 
(e.g. HIJING, apart form jet quenching). Also interestingly, the amount of Et released 
in the classical field model is twice the initial Et obtained in fl3| , so obviously a very 
strong energy loss would be needed in order to fall near the measurements of Et shown 
by PHENIX in this conference. 



3. Elliptic flow 

In non-central collisions, the anisotropic flow has been suggested clS cL SI mature of 
formation of pressure and as a probe of the early stages of the produced system [41]. 
Pressure can only be formed if the produced particles interact with each other. In non- 
central collisions the spatial azimuthal asymmetry of the production zone is transmitted 



s 



to the pressure gradients which drive the transverse flow. In non-central collisions, the 
initially generated flow is asymmetric, 7^ (Vy), leading eventually to an asymmetry in 
the azimuthal particle distributions. The harmonic coefficients of the Fourier expansion 
of the particle spectra are a quantitative measure of this asymmetry. In particular, the 
elliptic flow is characterized by the coefficient V2, defined as 
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Clear elliptic flow signal, up to 6% for v% and up to 15 % for Viipr) i n minimum bias 
events in Au-Au at y/s = 130 AGeV, has now been measured at RHIC by STAR [42 



(shown in Fig. [?[). An attempt to briefly summarize the situation with the semiclassical 
cascade models is presented in Fig. 6, where v% is shown from RQMDv2.4, UrQMD and 
ZPC at y/s = 200 AGeV (the measured v 2 can be expected to be close to that at y/s = 200 
AGeV [|7|): Partonic cascades can in principle generate large enough t> 2 but very large 
rescattering cross sections are needed As the mean-free paths of particles then 

become very short, this speaks in favour of hydrodynamics. Hadronic cascades predict 
the linear rise of V2{pr) at px ~ 500 MeV but they level off too soon after that [fl"4j| . These 
models do not seem to produce sufficient pr-averaged V2 either [0,0, which speaks in 
favour of an early formation of collectivity and pressure in the partonic system. 
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Figure 6. Elliptic flow for Au-Au at y/s = 200 AGeV from 
cascade models. From upper left: v% vs. b from RQMDv2.4 
5|, V2 vs. b from UrQMD f44j| , and V2(pt) vs. pt from 
UrQMD vs. px [44]. Lower panel: V2 vs. time from the 
ZPC parton cascade |f4l 



The observed strong elliptic flow confirms the prediction based on hydrodynamics [fl"7 



Detailed comparisons have also been done, as shown in Fig. [7|. The V2 obtained in the 
"conventional" hydrodynamic approach ||48|| , where the initial conditions are constrained 
by the observed final state multiplicity and where the Cooper-Frye procedure EB[ is 



applied for the decoupling, agrees with the data quite well. Similar conclusion is obtained 
in the Hydro-to-Hadrons approach H49 |, where the RQMDv2.4 hadron cascade is switched 
on below T = 160 MeV in order to treat the chemical and kinetic freeze-out in more 
detail. At very large impact parameters, where application of a purely hydrodynamical 
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system may be questioned in any case, these approaches seem to produce too much t> 2 . 
Also t>2<J>r) m the conventional hydrodynamic approach agrees with the STAR data very 
well up to 2 GeV or so. The deviation at high px is not surprising, since in a finite 
dynamic system with a finite lifetime the large-pr tails will not be adequately described 
by hydrodynamics. Another interesting observation is that the amount of Vi is not very 
sensitive to the equation of state (EoS) used. Therefore, the single particle spectra are 
needed to constrain the remaining uncertainties such as the EoS. 
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Figure 7. From left: V2 and v^ipr) from hydrodynamics 
right: vi from hydrodynamics followed by RQMDv2.4 [ 49 
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In conclusion, the first results from RHIC have again demonstrated the importance of 
a close interplay between theory and experiments: predictions for observables are difficult 
to obtain from truly first principles. Indispensable constraints for the models are obtained 
from certain sets of observables, such as the ones considered here. Especially useful will 
be the studies of the systematics in y/s, A and b. When the production and evolution 
stage of the strongly interacting system are theoretically under control, the signals of the 
QGP can be predicted and found. 
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